Chapter 2 


Quark-Gluon Plasma and the Early 
Universe 


There is now considerable evidence that the universe began as a fireball, the so called “Big-Bang”, 
with extremely high temperature and high energy density. At early enough times, the temperature 
was certainly high enough (T > 100 GeV) that all the known particles (including quarks, leptons, 
gluons, photons, Higgs bosons, W and Z) were extremely relativistic. Even the “strongly inter- 
acting” particles, quarks and gluons, would interact fairly weakly due to asymptotic freedom and 
perturbation theory should be sufficient to describe them. Thus this was a system of hot, weakly 
interacting color-charged particles, a quark-gluon plasma (QGP), in equilibrium with the other 
species. 

Due to asymptotic freedom, at sufficiently high temperature the quark-gluon plasma can be 
well-described using statistical mechanics as a free relativistic parton gas. In this Chapter, we 
explore the physics of QGP, perhaps the simplest system of strong-interaction particles that exists 
in the context of QCD. As the universe cooled during the subsequent expansion phase, the quarks, 
antiquarks, and gluons combined to form hadrons resulting in the baryonic matter that we observe 
today. The transition from quarks and gluons to baryons is a fascinating subject that has been 
difficult to address quantitatively. However, we will discuss this transition by considering the basic 
physics issues without treating the quantitative details. At present there is a substantial effort in 
theoretical physics to address this transition by using high-level computational methods known as 
lattice gauge theory. This subject is somewhat technical and we will discuss it only very briefly. 
However, the general features that have emerged from lattice studies to date are rather robust and 
can be discussed in some detail. 

The relatively cold matter that presently comprises everything around us is actually a residue 
of the annihilation of matter and anti-matter in the early universe. The origin of the matter- 
antimatter asymmetry which is critical for generating the small amount of residual matter is still 
a major subject of study, and we discuss this topic at the end of this Chapter. 

Another major thrust associated with the transition between the QGP and baryonic matter is 
the experimental program underway to study observable phenomena associated with the dynamics 
of this interface. This experimental program involves the collision of relativistic heavy ions that 
should produce (relatively) small drops of QGP. Large particle detector systems then enable studies 
of the products of these collisions, which can (in principle) yield information on the transition to 
the baryonic phase and the QGP itself. The program of experiments and the present state of the 
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experimental data will be discussed in Chapter XX. 


2.1 Thermodynamics of A Hot Relativistic Gas 


At very high-temperature such that the particles have energy much larger than their rest mass, 
we may describe them using relativistic kinematics and ignore their masses. Thus these energetic 
weakly interacting particles form a system that is, to an excellent approximation, a hot relativistic 
free gas. Since particles and antiparticles can be created and annihilated easily in such an environ- 
ment, their densities are much higher than their differences. Therefore the chemical potential u can 
be neglected. The number densities of the partons (species i) are then described by the quantum 


distribution functions B 
= Di 1 
n= PETT oy 
where 8 = 1/kpT and the — sign is for bosons and the + is for fermions. For relativistic particles, 
pi = E;. For Eib < 1, the exponential factor is small and there is a large difference between 
fermions and bosons. For £;3 > 1 the +1 becomes increasingly unimportant, and the distributions 
become similar. Integrating over the phase space, one finds, 


eres C(3)/T T? (boson) (2.2) 


(3/4)C(3)/n2T? (fermion) 


where ¢(3) = 1.20206... is a Riemann zeta function. The T?-dependence follows simply from 
dimensional analysis (the Boltzmann constant kg can also be taken to be 1). 

The energy density for a free gas can be computed from the same quantum distribution func- 
tions: 
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where the fermion energy density is 7/8 of that of boson. 
These expressions are valid for each spin/flavor/charge/color state of each particle. For a system 
of fermions and bosons, we need to include separate degeneracy factors for the various particles: 
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where gą = (9 + Zg) with g% and gy are the degeneracy factors for bosons and fermions, respec- 
tively. Each of these degeneracy factors counts the total number of degrees of freedom, summed 
over the spins, flavors, charge (particle-antiparticle) and colors of particles. When some species are 
thermally decoupled from others due to the absence of interactions (such as neutrinos at present 
epoch), they no longer contribute to the degeneracy factor. For example, at temperature above 100 
GeV, all particles of the standard model are present. At lower temperatures, the W and Z bosons, 
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top, bottom, and charm quarks freeze out and g, decreases. Therefore g, is generally a decreasing 
function of temperature. 

We can now calculate the contribution to the energy density from the quark-gluon plasma as a 
relativistic free parton gas. For a gluon, there are 2 helicity states and 8 choices of color so we have 
a total degeneracy of gẹ = 16. For each quark flavor, there are 3 colors, 2 spin states, and 2 charge 
states (corresponding to quarks and antiquarks). At temperatures below kgT ~ 1 GeV, there are 
3 active flavors (up, down and strange) so we expect the fermion degeneracy to be a large number 
like gf ~ 36 in this case. Thus we expect for the QGP: 


2 
eqcp © 47.5 (kT)! l (2.5) 


With two quark flavors, the prefactor is g, = 37. (For reference, if one takes into account all 
standard model particles, g, = 106.75.) 

The pressure of the free gas can be calculated just like the case of black-body radiation. For 
relativistic species, 


PH =e; (2.6) 
which is the equation of state. 
To calculate the entropy of the relativistic gas, we consider the thermodynamics relation, dE = 


TdS — pdV. At constant volume we would have just dE = TdS, or de = Tds where e (s) is the 
energy (entropy) per unit volume. Since e x T4, we can easily find that 


s= -=. (2.7) 
For an isolated system of relativistic particles, we expect the total entropy to be conserved. 
Now using Eq. 2.4 one can easily see that 
s x g (TT? , (2.8) 


where gą(T) counts the number of active (i.e., non-frozen) degrees of freedom in equilibrium. The 
total entropy of the active species is given by 


S x sR? «x g.(T)T?R? , (2.9) 


which is conserved in adiabatic processes. 


2.2 The Early Partonic Universe 


It has been established, since Hubble’s first discovery in the 1920’s, that the universe has been 
expanding for about ~ 10 billion years. The universe as we know it began as a “big bang” 
where it was much smaller and hotter, and then evolved by expansion and cooling. Our present 
understanding of the laws of physics allows us to talk about the earliest moment at the so-called 
Planck time tp ~ 10743 when the temperature of the universe is at the Planck scale T ~ My) 
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where Gy is Newton’s gravitational constant, and A and c are set to 1 unless otherwise specified. 
However, at this scale, the gravitational interaction is strong, the classical concept of space-time 
might break down. At times later than the Planck epoch when the universe has cooled below Mj), 
space-time may be described by a classical metric tensor g“”, and the laws of physics as we know 
them should be applicable. 

Since the observed universe is homogeneous and isotropic to a great degree, its expansion can 
be described by the Robertson-Walker space-time metric, 


dr? 
1 — kr? 


d3? = gu dz”dz” = dt? — R° (t) + r°(d8? + sin? 0de”) |, (2.12) 


which describes a maximally symmetric 3D space, where R(t) is a scale parameter describing the 
expansion and k is a curvature parameter with k = +1,—1,0, corresponding to closed, open and 
flat universe, respectively. 

The expansion of the universe after the Planck time is described by Einstein’s equation of 
general relativity, which equates the curvature tensor of the space-time to the energy-momentum 
tensor T””. The energy-momentum density comes from both matter and radiation and the vacuum 
Ag” contribution, the infamous “cosmological constant” of Einstein. If the matter expands as 
ideal gas, the energy-momentum density is 


Tuv = —P9uwv T (p + p)u”u” , (2.13) 


where p is the pressure and p is the energy density, and u” = (1,0,0,0) defines the cosmological 
comoving frame. The resulting dynamical equation for the scale parameter is 
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which is called the Friedmann (or Friedmann-Lemaitre equation). R/ R = H is the expansion 
rate (Hubble constant). Another equation needed for studying the expansion comes from energy- 
momentum conservation, 

p=—3H(p+p) . (2.15) 


Together with the equation of state p = p(p), the above equations can be solved to yield the 
evolution of p as a function of the scale parameter. There is now strong experimental evidence that 
we are living in a universe with k = 0 and A has been negligibly small until recently. Hence, we 
will focus below on a simplified Friedmann equation for the early universe without the second and 
third terms on the right-hand side in 2.14. 

When the temperature was lower than the Planck scale, the universe was an expanding gas 
of relativistic particles. These particles include quarks and leptons, the gauge bosons such as 
photons, gluons, and W and Z bosons, and perhaps more exotic particles like the supersymmetric 
partners of the standard model particles, heavy right-handed neutrinos, gauge bosons related to 
grand unification theories, etc. As the temperature cooled below the masses of certain particles 
(such as the W and Z bosons) they “freeze out” and decay, i.e., they are not longer created by 
inverse reactions of their decay products due to the lower temperature. Some of these particles 
with a short life time had disappeared long ago, and some with a long life time may still be with 
us today in the form of dark matter. 
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Thus we expect that when the temperature drops below the electroweak scale (T < 100 GeV) 
the early universe will be a hot gas of the standard model particles: quarks, leptons, gluons and 
photons. Since the system is dominated by the strongly interacting degrees of freedom, quarks 
and gluons (i.e., partons), it is a good approximation to regard it as a system of quark-gluon 
plasma. Because of asymptotic freedom, the interaction between quarks and gluons are fairly weak 
at high-temperature, and it shall be a good approximation to describe the plasma in terms of a 
non-interacting parton gas. 

During this phase of the universe, the energy density p is dominated by these relativistic partons 
and decreases as the universe expands. The evolution of p during this time is governed by the fact 
that we have a gas of relativistic partons. The volume of any piece of the universe increases like 
R?, but the energy in every mode decreases as R~! (as the wavelength of the mode expands with 
the universe). Thus we expect 

pxR*, (2.16) 


and Eq. 2.14 then yields 
Re Ro, (2.17) 


which has the solution R ~ vt. That is, the size of the universe increases as the square root of 
time. The energy density then decreases as p ~ t~?. 

If we assume that the number of effective degrees of freedom, gą is constant during the early 
evolution of the radiation-dominated universe then the radiation energy density (pọ x Tt, as in 
Eq. 2.4) with its variation as R (Eq. 2.16), we find that the temperature varies inversely as the 
radius parameter T x R`! and therefore T x t~!/?. Note that according to Eq. 2.7 this also implies 
that the total entropy of the universe is conserved. We then obtain the following relation for the 


temperature as a function of time: 
AMp, 
T(t) > lesz ire (2.18) 


(g4) 2T2 


If we invert this relation to yield 


te (2.19) 


we can construct the timeline for the temperature of the early universe from 10743 sec. through 
about 10° yr. when the radiation dominated phase ends. 


molecules, 
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Figure 2.1: History of the universe for temperatures less than kgT ~ 100 GeV. 
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We have assumed that g, is constant in obtaining these results. However, we do need to consider 
the fact that as the temperature drops some particles freeze out, and so g,(T') then changes. This 
will modify the expressions 2.18 and 2.18. However, the basic behavior of the expanding universe 
is qualitatively described by these relations, especially noting that in Eq. 2.18 the dependence of 


the temperature on gx is very mild (T x gy a D 


2.3 The Quark-Gluon Plasma in Perturbative QCD 


Until this point we have been treating the quarks and gluons in the QGP as free particles without 
interactions. Of course, in a high-temperature QGP we expect QCD perturbative theory to be ap- 
plicable due to asymptotic freedom. One important additional consequence is that chiral symmetry 
is now a good symmetry, and the chiral condensate must vanish in the plasma 


(QGP|IYY|QGP) = 0, (2.20) 


where strictly-speaking the QGP “state” is actually the thermal average over the excited states of 
the QCD vacuum when the baryon number density is ignored. 

Another important feature of the QGP is color deconfinement. In QCD perturbation theory, 
the quarks and gluons are free particles that can be described by plane waves. Asymptotic freedom 
will guarantee that high momentum transfer interactions are weak. Small-momentum transfer 
scattering involves long distance interactions which are screened by the plasma (although this is 
only strictly true in the color electric sector). As such, the charged quarks and gluons can move 
freely inside the plasma without being confined to a local region. This remarkable property is 
radically different from the low energy limit of QCD where all charges are permanently confined to 
the interior of hadronsa scale about 1 fermi. 

Consider a color charge in midst of a color-neutral plasma. The other particles in the plasma will 
act to screen it, and as a consequence the interaction between color charges is damped exponentially. 
To calculate the screening length, one can start from a color charge and calculate its induced color 
fields. The result is a correlation function of gluon fields. This function can be calculated in 
perturbation theory at high-temperature, and the result for the screening mass is 


ee ee (2.21) 


to leading order in the strong coupling expansion. The so-called Debye screening length is simply 
1/mp or 1/gT, which is very short at high-temperature. When the color charges are screened in a 
plasma, it has a finite energy and therefore in this sense, the color charges are now liberated. 

Unfortunately, the magnetic interaction is only weakly screened; it has a screening mass of order 
g°T. Absence of the magnetic screening means that the magnetic sector of the QCD remains non- 
perturbative even at high-temperature. Fortunately, at high-temperature this non-perturbative 
part contributes to physical observables only at higher-order in QCD coupling, so the free gas 
behavior is dominant. 

Another important feature of the plasma is the plasma frequency. In a QED plasma, light 
cannot propagate below the plasma frequency, wp) = (ne? /m)}*/ 2 but will be reflected from the 
surface, like in a silver-plated mirror. The physics of the QGP is similar: gluons (plasmon) cannot 
propagate as a free field in the plasma if its energy is too low. In fact, the gluons acquire an effective 
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mass which is effectively the plasma frequency. Perturbative calculations confirms this behavior, 
and to leading order in perturbation theory the plasma frequency is 


T EJN. + N;/2(gT) . (2.22) 


where Ne = 3 is the number of color and Ny is the number of fermion flavor. The transverse- 
polarized gluon modes acquire the same mass. 

The plasmon and transverse gluon modes are damped in the plasma. One can calculate the 
damping rate using the so-called hard-thermal loop method in pQCD and the result is gauge- 
invariant: y = ag*NcT/(247), where numerically a is found to be a = 6.63538. 

The results we discussed above are the basic leading-order predictions of pQCD. Higher-order 
contributions can and have been calculated in the literature. Unfortunately perturbative expansions 
for thermodynamical properties of the plasma converge very slowly. As such, the free plasma picture 
works only at extremely high-temperature. Even at the temperatures corresponding to 100 GeV 
the perturbative expansion must be reorganized significantly to get a sensible prediction. We will 
come back to this point later. 


2.4 Transition to the Low-Temperature Phase: Physical Argu- 
ments 


As we have discussed in the previous chapter, the zero temperature ground state of QCD is strikingly 
different from the high-temperature QGP: color charges are confined to the interior of individual 
hadrons and chiral symmetry is broken spontaneously. Therefore, as the plasma cools in the 
universe, some rapid changes in thermodynamic observables must occur from the high-temperature 
QGP phase to the low-temperature confining and chiral-symmetry breaking phase, where the quarks 
and gluons combine to form colorless states of hadronic matter. 

It is possible to estimate the transition temperature by comparing the QGP gas pressure with 
that of hadronic gas. The lightest hadrons are pions, and for T < 1 GeV (note that in the following 
we often use units where kg = 1 and T has units of energy), we might expect a gas of relativistic 
pions. This is a system with only 3 degrees of freedom, g = 3, so the energy density and pressure 
of the system 

= 3T a P, = —~—T* (2.23) 
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This, however, is not the full story. Pions are collective excitations of the non-perturbative QCD 
vacuum. This true ground state of the QCD vacuum has a lower-energy —B than the perturbative 
QCD vacuum. (In the MIT bag model of hadrons, this energy is the origin of the quark confine- 
ment.) Lorentz invariance requires that the energy-momentum density is of form TH” = Bg"”. Thus 
the non-perturbative QCD vacuum has a positive pressure as well. Therefore, the total pressure of 


the hadronic phase is 
Pow = B+ 3a 4 ; (2.24) 
90 
On the other hand, from the previous sections, the pressure of the QGP phase with 2 quark flavors 


is, Pagp = 37n°T*/90. Equating the two pressures, we find the transition temperature, 


T, = (45B/171?)\/4 ~ 180MeV , (2.25) 
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where we have used the MIT bag constant B = 200 MeV as determined by fits to the masses of 
physical hadrons. 

The energy difference (latent heat) between the two phases at the transition temperature is 

34r? 
Ap= 2 4B, (2.26) 
30 

which is on the order of 2 GeV/fm?. 

Another estimate of the transition temperature comes from considering chiral symmetry. At 
finite but small temperature, the pion gas will dilute the chiral condensate in the zero-temperature 
vacuum. The quark condensate can be calculated as a response of the system’s free energy to the 


quark mass, 
1 OF 


(hr = Ny ðm , (2.27) 


where Ny is the number of light quark flavors. The free-energy of the pion gas is 


F = (N?- yr f In(1 — e~Bs/T) , (2.28) 


Thus the pion condensate has the following low-temperature expansion, 


N2-1T? 3 
BN, 4f2 0 "|? 


(yir = po Í (2.29) 


where the ellipse indicates higher-order terms in the expansion. If one just keeps the first two 
terms, the chiral condensate vanishes when 


Te = 2fny/3Ny¢/(N? — 1) = 200 MeV , (2.30) 


which is consistent with the other estimate. 

Clearly, the QCD system is strongly interacting around Te. On the other hand, the above 
estimates relied on calculations which are valid at temperatures much higher than Te. To say 
something rigorous about what happened around Tą, one must resort to lattice QCD, a numerical 
approach to solve QCD through computer simulation. 


2.5 A Brief Tour in Lattice QCD Thermodynamics 


At lower temperatures where the coupling constant is larger one must employ non-perturbative 
methods of calculation. The only known method for solving QCD non-perturbatively is on a space- 
“time” lattice. Here we present a simple introduction to this method without getting involved in 
too many technical details. 

Consider a QCD system with temperature T = 1/3 and baryon number zero (this is true to 
a good approximation in the early universe, as we will discuss later in this chapter). The most 
important quantity is the partition function, 


Z = Trlexp(—GA)] , (2.31) 
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where H is the QCD hamiltonian and trace is over all physical states in Hilbert space. H is a 
function of 3D quark and gluon fields y(7) and A“(#), respectively. One can introduce a fourth 
coordinate x4 (imaginary time) and the 4D field 


OL, z4) = efta g(@)e— Hea ; (2.32) 


where @ collectively labels all QCD fields and x4 runs between the values 0 and 8. Then the 
partition function can be written as a path integral 


Z = Tr | Did|exp(-Sl¢)) . (2.33) 


where D[|¢] is a functional integration measure and S[A] is the QCD action defined by the lagrangian 
in Euclidean time, 


S{¢] = i des | PELE. : (2.34) 


Therefore, the partition function is now reduced to a functional integral. 
Because of the trace in Z, the gluon potentials A“ obey the period boundary condition in x4, 


A" (x4, 2) = A" (x4 + 8,2) , (2.35) 
whereas the quark fields obey the anti-periodic condition 
o" (x4, %) = —b" (x4 + 8,2) . (2.36) 


The whole thermodynamics formulation is then invariant under SU(3) gauge transformation with 
U(a4,£) satisfying the periodic boundary condition. 

The path integral formulation makes explicit that QCD is a quantum field theory with an 
infinite number of quark and gluon degrees of freedom. To solve it approximately, one first limits 
the system to a 3D box of dimensions Lg, Ly, Lz, and then replaces the continuous space- “time” 
variables by a discrete 4D lattice. The quark and gluon degrees of freedom now live on the lattice 
sites and the bonds in-between the sites, respectively. The number of these d.o.f. is now finite but 
large (millions to billions in actual simulations). 

The integrations over quark fields are of Gaussian type and can be done analytically. Those over 
the gluon fields can be evaluated using the Monte Carlo sampling method. In this approach, classical 
gluon configurations are generated on the lattice with the probability distribution corresponding 
to the Boltzmann factor e79. The actual physical observables are calculated with hundreds and 
thousands of these “typical” distributions. As an example, the Debye screening mass can be 
measured from the space correlation of the gluon fields at equal time. One of the important tools 
of the lattice calculation is that one can change the parameters of the theory and study how 
observables change in a world different from the real one. 


2.6 Quark Masses and The Nature of the Transition 


The nature of the QCD transition from high to low temperature depends on strongly on light-quark 
masses. For simplicity, let us ignore heavy quarks (c, b, t) and concentrate on up, down and strange 
light quarks. 
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One powerful theoretical approach to discuss phase transitions is Laudau-Ginsburg theory. In 
this framework, one first identifies the order parameter of a transition, which is the fundamental 
observable which drives the phase transition. An example of order parameter in water-steam 
transition is the density. One then constructs an effective lagrangian of the order parameter which 
governs the dynamics of the transition. Here the fundamental symmetry constraints are important 
to determine the form of the effective lagrangian, from which one can often say a lot about the 
transition without going into the details of the specific dynamics. 

If all three quarks are massless, the low-energy order parameter for chiral transition can be 
chosen to be ©, a 3 complex matrix which transforms under chiral SU(3), x SU(3)p as 


o> UpdUE, , (2.37) 


One can then construct an effective lagrangian to describe the dynamics of X. One of the term 
that one can write down is Det» which is trilinear in components of the order parameter. As a 
consequence, the chiral phase transition described by the theory is in general a first-order phase 
transition. 

When the strange quark is infinitely heavy, and up and down quarks remain massless, the 
low-energy order parameter can be taken to be a 2 x 2 unitary matrix. The phase transition in 
this system is similar to that of an O(4) magnet, a magnetic with four-independent magnetization 
direction. The Landau-Ginzburg theory for this system leads to a second-order phase transition. 
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Figure 2.2: Order of phase transition as a function of light quark masses (mu = ma) vs. the strange 
quark mass(m,). 


Therefore, if one varies the mass of the strange quark from large to small, the second-order 
phase transition must end at a tri-critical point, beyond which the transition becomes first order. 
The tricritical point is characterized by vanishing of the coefficients of the first two terms (quadratic 
and quartic) in the potential for the order parameter. 

The above consideration is not yet sufficient because the up and down strange quark masses do 
not vanish in the real world. 
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If up, down and strange quark masses are taken to infinity, we are left with a theory with just 
pure glue fields. This system at finite temperature has an interesting symmetry called color Z3 
center symmetry. An effective Landau-Ginzburg theory can be constructed to describe the physics 
of Z3 symmetry breaking at high-temperature phase. The symmetry argument suggests that the 
transition is of first order. The strength of the first-order phase transition lowers as the quark mass 
becomes lighter, and finally the first-order phase transition region is enclosed by a second order 
phase transition line. 

A diagram detailing the nature of QCD phase transitions as a function of quark masses is shown 
in Fig. 2.2. In the limit of 3 very light quarks, we expect a first order transition. The first order 
transition region is enclosed by a second order transition line which goes through the tri-critical 
point. In between the two second-order phase transition lines, one has a broad region of rapid 
cross-over. The exact locations of these second-order transition lines are not known. Therefore, for 
a realistic physical situation, where the up and down quark masses are on the order of a few MeV 
and strange quark mass is about 100 MeV, the transition can either be of a weak first-order or a 
rapid crossover (2nd order), as shown by two dots in Fig 2.2. 


2.7 Physics of the QCD Transition on Lattice 


Lattice simulations of QCD thermodynamics have made significant progress in the last decade, 
due to both rapid rise in computational power and implementation of better algorithms. Some of 
challenges in achieving a complete realistic simulation include finite lattice size effect, discretization 
errors, implementing dynamical quarks, and simulations at small quark masses. 

From simulations of QCD on a lattice, a transition in thermodynamic observables is clearly seen 
at a fairly well-defined temperature of about To ~ 150 MeV. The energy density undergoes a rapid 
change near a critical temperature Tc, enhanced by almost an order of magnitude, as indicated in 
Figure 2.3. Beyond Tç, the energy density is fairly flat as a function of temperature but slightly 
below that predicted by the free gas model. This rapid change is an indication that the fundamental 
degrees of freedom are different above and below Te. 

The transition is less dramatic in the equation of state, i.e., the pressure of the system as a 
function of temperature. At low-temperature, the pressure is very small. As temperature increases, 
the pressure builds up gradually over a large range of T, from Te to 2T,. When the pressure curve 
flattens out at high-T, it again undershoots the result of the free gas model. The equation of state 
is also more sensitive to different quark mass scenarios. These calculations are improved constantly 
with smaller quark masses and lattice spacing. Improved calculations show that the transition at 
physical quark masses is a rapid cross over. 

An interesting property of this transition emerges from the consideration of chiral symmetry in 
lattice simulations. Recall that, in the absence of quark masses, the QCD Lagrangian is chirally 
symmetric, i.e., nature is invariant under separate flavor rotations of right and left-handed quarks. 
This symmetry is evident in the QGP phase at T > Tce in the lattice simulations. This is quite 
analogous to the symmetric demagnetized state of a ferromagnet above the Curie temperature. 
The ferromagnetic property is related to the magnetization M which vanishes above the Curie 
temperature. The lattice QCD calculations can provide a measurement of the relevant “order 
parameter”, the scalar quark density (Yy), as a function of temperature. This order parameter 
provides a measure of the effective mass of a quark in the medium. The result is that at T > To 
this effective quark mass becomes small, approaching zero as T — oo, as one would expect in 
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Figure 2.3: The transition from mesonic matter to the QGP phase as suggested by lattice simu- 
lations of QCD. The simulations are carried out with 2 or 3 light flavors or 2 light and 1 heavy 
flavor. The expected limits as given by Eq. 2.5 are shown by arrows on the right side of the figure. 
(Note that in this figure T actually represents kg times temperature.) 


chirally symmetric QCD. As shown in Figure 2.4, below Tg there is a sharp increase in (ypy) 
corresponding to the quarks developing a constituent mass of ~ 300 MeV. This heavy consituent 
quark is the basis of the quark model of hadrons to be discussed in Chapter XX. This transition 
to the broken symmetry phase is again analogous to the ferromagnet, which spontaneously breaks 
rotational symmetry in developing a finite M below the Curie temperature. Shown in the same 
figure is the susceptibility of the condensate which shows a peak at the transition temperature. 
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Figure 2.4: The scalar quark density (py) measured in lattice QCD as a function of temperature. 
At T < To chiral symmetry is broken, whereas the chirally symmetric limit is realized at higher 
temperatures T > To. 


The potential between two quarks is shown for various temperatures in Figure 2.5. One can 
see that at low temperatures T < Tc, the potential continues to rise at larger distances, consistent 
with the expectation that the quarks will be confined as colorless hadrons. At higher temperatures, 
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T > Tg, one finds that the potential energy at large distances saturates, and it is possible for the 
quarks to propagate as a free particle. This confirms the picture that above Tc one has a deconfined 
plasma. 
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Figure 2.5: The potential energy between two quarks as a function of the separation distance as 
computed by lattice QCD. The calculation is performed for various temperatures showing that 
confinement is present at low temperatures T < Te but not above To. 


2.8 Evolution of the Universe in Hadronic Phase 


The QGP dominated universe will expand and cool until we reach the critical temperature, Tc, to 
produce a gas of pions along with a few baryons and antibaryons. The number of baryons is clearly 
suppressed by the Boltzmann factor 

PweMa/T | (2.38) 


which is small but still significant. But the universe we observe is clearly not a pion gas, nor the 
decayed remnants which would be electrons, positrons and finally just photons. As the universe 
cools, the pion densities diminish exponentially as a function of temperature. This decrease is 
achieved by the annihilation 

nET 3y+)7, (2.39) 


whereas the inverse reaction is difficult because the thermal photons do not have enough energy 
to produce pions. Similar process happens for baryons and antibaryons, while they continue to 
annihilate, they cannot be reproduced. If the annihilation rate is rapid enough, the density of 
baryons follows the thermal Boltzmann distribution. Of course, if the universe had no net baryon 
number, all baryons would disappear eventually. However, the existence of matter around us shows 
that the baryon number of the universe is not zero. Therefore, as temperature cools, all anti-baryons 
are annihilated, all pions are either annihilated and/or decayed, but there is a small tiny baryon 
density in the form of proton and neutron survives. (The origin of these baryons is the subject of 
the next section of this chapter. Indeed, these protons and neutrons are the main characters of the 
remainder of this book.) 

What happened to other components of the universe? As temperature lowers, we are left with 
just leptons and photons. At temperatures below the electron mass, pairs of electrons and positrons 
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are no longer created, so they freeze out and annihilate to produce more photons. We enter a phase 
of the universe that is dominated by photons and neutrinos: a black-body universe. Of course, 
due the charge neutrality, a small fraction of electron residue is also present. The energy density 
will continue to decrease as t~? until much later when, it turns out, the small fraction of baryonic 
matter becomes a significant factor in the energy density. 

We still observe the remnant black-body radiation associated with the early universe. (The 
neutrinos have not been observed, but are discussed in Chapter 8.) There is presently a uniform 
distribution of cosmic microwave radiation with a characteristic temperature of (2.725 + 0.001)° K 
and a number density of 410.4+0.5/cm%. The present baryon to photon ratio is tiny, about 6x 107!°, 
but the photons have energies of ~ 2.5 x 1074 eV and the baryons have energy Mp ~ 1 GeV. Thus 
the energy density of the baryons is about 2400 times the energy density of the photons. As the 
universe cooled from the “radiation dominated” era to the present, the photons were red-shifted 
due to the expansion. The photons decoupled from the protons and electrons when the latter 
combined into neutral hydrogen. Since then (temperature T ~ 0.3 eV) the photons have been 
red-shifted by about 10*, while the protons remained at M, ~ 1 GeV. Thus we now find ourselves 
in an era where baryonic matter dominates over the cosmic photons. (In addition, we have dark 
non-baryonic matter and dark energy which remain to be understood - but that is a different story 
and we will confine our attention here to the baryons and photons.) 


2.9 The Origin of Baryonic Matter 


Let’s now return to the source of the baryonic matter. While the present baryon to photon ratio 
is indeed tiny, that is the stuff of which we are made. Since the early universe should consist of 
equal numbers of quarks and antiquarks, eventually they should all annihilate to produce only 
photons in the end. The observed number of baryons is much greater than can be expected from 
a random fluctuation associated with the cooling of the enormous number of quarks and gluons in 
the initial plasma. The annihilation of baryons and antibaryons ceases at a rather low temperature 
of ~ 20 MeV. The number density of baryons and antibaryons at this temperature is given by the 
expression (for gg species of non-relativistic particles) 


MT 3/2 _ Mp 
nB = 9B (==) eT, (2.40) 
The corresponding number density of photons is 
2 
ia aa (2.41) 
so the ratio is approximately 
nB Mp 3/2 _ Mp 
pets | AG E T 2.42 
Ny oe ( T ) i (2.42) 
~ 107'% (@T = 20 MeV). (2.43) 


We would expect the chance excess of baryons over antibaryons to be a very small fraction of the 
total number of both. Therefore, the observed 6 x 1071? fraction of baryons relative to the number 
of photons is not possible to generate by just a statistical fluctuation. 
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So what is the mechanism that generates the small excess of matter over antimatter? The 
modern view is that the observed excess of baryonic matter relative to antimatter is likely due 
to the properties of particles and their interactions during the expansion and cooling of the early 
universe. This can occur if three criteria, initially discussed by Sakharov, are satisfied: 


1. baryon number (B) is violated so that baryons (and antibaryons) can be created, 
2. CP is violated so that the rates of baryon and antibaryon production can be different, 
3. thermal equilibrium is broken so that forward and reverse reactions become unbalanced. 


Studying baryogenesis has been a very active area of theoretical physics for the last few decades. 
In the standard model, the baryon number can be violated through non-perturbative process via 
the so-called anomaly. However, to generate enough baryons, one must have a strong first order 
phase transition which is hard to achieve with the heavy Higgs mass. Moreover, the standard model 
CP violation through the CMK matrix is too small to generate enough CP asymmetry. Therefore, 
it appears that baryogensis requires physics beyond the standard model. One of the interesting 
direction is what happens in a supersymmetric extension of the standard model. 

In grand unified models, baryon number is in general violated. One consequence is that protons 
would be unstable and decay, typically by a process like p > e+7°. Searches for such decays limit 
the proton lifetime to Tp > 103? years. Consider a heavy particle of mass My. Fermi’s golden rule 
can be used to calculate the decay rate as 


r = |(f| Axe) (2.44) 
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The matrix element will be proportional to ax /M? where ax is the squared coupling constant 
and the M iy is the propagator for the heavy particle exchange. Since the only other energy scale 
in the problem is the proton mass M, we can write by dimensional analysis 


2 
ee (=) . (2.45) 


Using the experimental limit for Trp and assuming ax ~ 0.1 one obtains a rough limit for the 
mass of the heavy particle responsible for baryon number violating interactions Mx > 10!° GeV. 
In a typical baryogenesis scenario, one might generate long-lived heavy particles X in the very 
early universe when T > Mx. At T < Mx the X particles freeze out and then much later, when 
the universe is cooler, they decay into baryons and antibaryons. The presence of CP violation 
in the decays enables a production of an excess of baryons over antibaryons. However, it is now 
generally believed that this excess gets washed out at lower temperature by the so-called electroweak 
sphaleron process. 

In more recent years, massive right-handed neutrinos have been considered as good candidates 
to lead to baryogensis. These particles could have masses slightly smaller than the grand unifi- 
cation scale, and decay asymmetrically (due to CP violation in lepton sector) to produce a net 
lepton number. This lepton number can be converted in to baryon number through the sphaleron 
processes. This route of generating baryon asymmetry is called leptogenesis. 
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2.10 ***Appendix for Chapter 2**** 


2.10.1 Ginzburg-Landau Theory for Phase Transitions 


In a first-order phase transition, the free-energy is continuous and its first-order derivatives are 
not. In a second-order transition, the first-order derivative of the free energy is continuous and the 
second-order derivatives are not. 

Quite often, phase transitions involve a transition from a symmetrical phase to a less symmet- 
rical one, or vice versa, i.e., a symmetry-breaking transition. For instance, in the transition from 
a paramagnetic to a ferromagnetic system, the rotational symmetry is broken because a sponta- 
neous magnetization defines a unique direction in space. In transition from normal liquid 4He to 
superfluid liquid 4He, gauge symmetry is broken. Near the critical point in liquid-gas transition, 
the distinction between liquid and gas disappears above the critical point. 

Because of this rather distinct feature of the transitions involving symmetry-breaking, a new 
macroscopic parameter was introduced by Landau to describe the transition phenomenologically, 
and is called order parameter. The order parameters take zero in a symmetrical phase and non- 
zero in the unsymmetrical (or less symmetrical) one. The order parameters may be a scalar, vector 
or tensor, a complex number, or some other quantity, depending on the symmetry of the transition 
involved. The order parameter changes continuously near the second phase transition point, so the 
volume or entropy do not change abruptly. For this reason, a second-order phase transition is also 
called a continuous transition. One important difference between the order parameter and other 
macroscopic variables such as pressure and temperature is that the values of the order parameter 
are determined by minimizing the thermodynamic potential of a system. 

Ginzburg and Landau found a general way to describe symmetry-breaking phase transitions 
in terms of a free energy functional involving order parameters. For simplicity, let us assume the 
order parameter is a vector 7 and construct a free energy which has a minimum at 77 = 0 above 
the transition point (T > Tc) and 7 4 0 below it. The free energy which is a scalar function of 
the order parameter, depending on the scalar-scalar product 77-77. Near the phase transition point 
where |7| is small, one can make the following Taylor expansion, 


S(T, if) = S(T) + a2(T) |i? + aa +> (1) 


If we choose a2(T) = ao(T — Te), then when T > Tc, ag > 0 and 7 = 0 is a local minimum of ®, 
when T < Tc, ag < 0 and 7 = 0 is a local maximum. This can be seen by plotting ® as a function 
of |7| near 7 = 0. This choice makes the order parameter behave in the way described above. To 
ensure 7 = 0 is also a global minimum for T > Tc, we take a4(T) > 0 at all T. 

If we neglect the high-order terms in (1), the potential and the order parameter of the system 
evolve with temperature in the following way. At T > To the potential is shown in Fig. xx: 7 = 0 
is the minimum and the system is in the symmetrical phase. At T < Tg, the potential is shown in 
Fig. 1b, and there are minima at |7| = const. with arbitrary phase and a local maximum at 7 = 0. 
These can be obtained from 


® 
OF = agi + das [P F=0 , (2) 
which gives 
3 - — Q2 
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The second equation tells us that the order parameter changes as (Te — T)? with 8 = 5 below Tg. 
G is one of the critical exponents that are introduced to characterize the singular behavior of an 
observable near the critical point. We will introduce more critical exponents below. 

Substituting (3) into (1), we find the free energy 


7 ®o(T) T > Te a 
(T) = 4 
(T) — a2 / fag T< Te 


Thus ® and its first derivative are continuous across To. However the second-order derivative which 
is related to specific heat, 


ao 
C = -T— 
aT? 65) 
is discontinuous. It is easy to show 
ag 
Clrer; = Clr-r4= TA (6) 


If one defines C œ |T — Te|7® where a is another critical exponent, then a = 0. 

When an external field h is applied to the system, the potential is added with a term řī- hV, 
where V is the volume of the system. Then 77 4 0 at any temperature, and the second-order phase 
transition disappears. Eq. (2) becomes 


2a + 4a4lq]? 7 = hV (7) 


Above To, |7j| is small, and we have 7 = AV / 2a. The susceptibility is 


On V 2 
x=] =5— ~ |T- To (8) 
Ohln=0 202 
where the critical exponent y is 1. Below To, 7 = hV/(—4a2) and again y ~ |T — To|~7 with 
/ 
yi = 1. 
At T = Tc we have the following relation between order parameter and the applied field, 
7 (ee ie a @) 
n= da , 


where the critical exponent ô is 3. 

The critical exponents in Landau theory are independent of the specific values of the parameters 
that appear in (1). They are completely universal. In other words, the phase transitions in very 
different physical systems exhibit the same singular behavior. A general theory of phase transition 
must capture this universality feature. 


2.11 Problem Set 


1. Calculate the number and energy densities of a relativistic fermion/boson gas. And show that 
in the standard model gą = 106.75. 
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2. Calculate the temperature as a function of time during radiation-dominated era of the early 
Universe, up to the normalization factor. 

3. Calculate the free-energy of a free pion gas, from which derive the dependence of the chiral 
condensate in temperature to first order in T?/ f2. 

4. Show that the thermodynamic function of QCD without quarks has a Z3 symmetry. Construct 
a Ginzburg-Landau theory for Z3 phase transition and show it is a first-order transition. 

5. Calculate the baryon energy density over the photon energy density in the universe. 


